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ABSTRACT 

The origin of the X-ray background, in particular at hard (2-10 keV) energies, has been a debated issue for more 
than 30 years. The Chandra deep fields provide the deepest look at the X-ray sky and are the best dataset to study 
the X-ray background. We searched the Chandra Deep Field South for X-ray sources with the aid of a dedicated 
wavelet-based algorithm. We are able to reconstruct the Log N-Log S source distribution in the soft (0.5-2 keV) 
and hard (2-10 keV) bands down to limiting fluxes of 2 x 10~ 17 erg s _1 cm~ 2 and 2 x 10 -16 erg s _1 cm~ 2 , 
respectively. These are a factor ~ 5 deeper than previous investigations. We find that the soft relation continues 
along the extrapolation from higher fluxes, almost completely accounting for the soft X-ray background. On the 
contrary, the hard distribution shows a flattening below ~2x 10~ 14 erg s _1 cmT 2 . Nevertheless, we can account 
for ,> 68% of the hard X-ray background, with the main uncertainty being the sky flux itself. 

Subject headings: diffuse radiation - surveys - cosmology: observations - X-rays: general 



1. INTRODUCTION 

In the soft (1-2 keV) band, ROSAT surveys resolved 70 - 
80% of the cosmic X-Ray Background (XRB) into discrete 
sources, at a flux level of 1 x 10~ 15 erg s _1 cm -2 (Hasingeret 
al. 1998). The great majority of sources brighter than 5 x 10 -15 
erg s" 1 cm -2 were optically identified with unobscured Ac- 
tive Galactic Nuclei (AGN; Schmidt et al. 1998; Lehmann et 
al. 2001). The origin of the background at higher energies is 
more controversial. The ASCA and BeppoSAX deep surveys 
achieved relatively shallow detection limits of (3 — 5) x 10 -14 
erg s _1 cm -2 , resolving ~ 25% of the XRB into point sources 
(Ogasaka et al. 1998; Cagnoni et al. 1998; della Ceca et al. 
2000; Giommi et al. 2000). Main contributors are thought to 

' be absorbed and unabsorbed AGN with a mixture of quasars 
and narrow emission-line galaxies as optical counterparts (e.g. 
Fiore et al. 1999; Akiyama et al. 2000). Recently, Chandra 
and XMM-Newton have dedicated long exposures to address 
this issue (Mushotzky et al. 2000; Giacconi et al. 2001a; Horn- 

■ schemeier et al. 2000, 2001; Brandt et al. 2001; Hasinger et 
al. 2001; Tozzi et al. 2001). The analysis of the first 300 ks of 
the Chandra Deep Field South (CDFS) allowed to reach limit- 
ing fluxes of - 10" 16 erg s" 1 cm~ 2 (0.5-2 keV) and - 10~ 15 
erg s _1 cm~ 2 (2-10 keV; Tozzi et al. 2001). Similar results 
were obtained on the Hubble Deep Field by Chandra observa- 
tions (Hornschemeier et al. 2000, 2001). These observations 
allowed the resolution of 81 - 95% of the soft (1-2 keV) XRB, 
confirming the results from the deep ROSAT HRI surveys in 
the Lockman hole (Hasinger et al. 1998). In the hard band, the 
absolute value of the background is poorly established and, de- 
pending on the adopted normalization, Tozzi et al. (2001) were 
able to account for 60 - 90% of the hard (2-10 keV) XRB. 

The increasing angular resolution and source crowding pro- 
vided by the current generation of X-ray instruments demands 
advanced detection algorithms. One of these is the wavelet 
transform (WT), a mathematical tool for multi-scale image 
analysis. The automatic detection and characterization of 
sources in X-ray images by means of WT-based algorithms 
were first carried out by Rosati et al. (1994, 1995). WT- 
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based techniques have been developed and applied also by other 
groups to the analysis of ROSAT images (e.g. Slezak et al. 
1994; Grebenev et al. 1995; Damiani et al. 1997; Vikhlinin et 
al. 1998; Lazzati et al. 1999). In view of the excellent results 
so far obtained, it is natural to extend these techniques to the 
analysis of the far richer Chandra images (Freeman et al. 2001; 
Harnden et al. 2001: Moretti et al. 2001). Here we analyzed 
the full (0.97 Ms) CDFS (Giacconi et al. 2001b; Rosati et al. 
2001) by means of an advanced detection algorithm originally 
developed for the analysis of the ROSAT HRI fields (Lazzati 
et al. 1999; Campana et al. 1999) based on the WT (Brera 
Multiscale Wavelet - BMW). At variance with other WT algo- 
rithm, our BMW characterizes detected sources by means of a 
decimation scheme (Lazzati et al. 1999). The application and 
testing of this algorithm on Chandra fields and on the CDFS in 
particular is extensively described in Moretti et al. (2001). 

2. CHANDRA OBSERVATIONS 

The CDFS consists of eleven observations for a total expo- 
sure time of 0.97 Ms. All exposures were taken with the Chan- 
dra X-ray Observatory (Weisskopf et al. 2000) Advanced CCD 
Imaging Spectrometer (ACIS-I) detector (Garmire et al. 2001). 
ACIS-I consists of four CCDs arranged in a 2 x 2 array for a 
field of view of 16.9' x 16.9'. The on-axis image quality is 
- 0.5" FWHM, increasing to - 3.0" FWHM at - 4.0' off- 
axis. The data were filtered to include only the standard event 
grades 0, 2, 3, 4 and 6 (see Chandra proposers' observatory 
guide). All hot pixels and bad columns were removed. We 
removed also flickering pixels with more than two events con- 
tiguous in time (time interval of 3.3 s) and residual cosmic ray 
events. 

The eleven observations were co-added with a pattern recog- 
nition routine to within 0.5" r.m.s. pointing accuracy (Moretti 
et al. 2001). We restricted our analysis to the ACIS-I area fully 
exposed in the eleven observations, selecting a circle with 8' 
radius, for a total coverage of ~ 0.056 deg 2 . In order to com- 
pare our results with previous studies we carried out the detec- 
tion in two selected energy bands: a 0.5-2 keV soft band and 
a 2-7 keV hard band (in the range 7-10 keV the effective area 
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steeply decreases while the background increases resulting in 
a lower signal to noise ratio for X-ray sources; fluxes in the 
hard 2-7 keV band are extrapolated to the classical 2-10 keV 
energy band). Time intervals during which the background rate 
is larger than 3 a over the quiescent level were removed in each 
band separately. This results in a net exposure time of 942 ks 
in the soft band and 936 ks in the hard band. The average back- 
ground over the fully exposed portion is 0.14 (0.19) counts per 
square arcsec in the soft (hard) band. Following Tozzi et al. 
(2001), the count-rate to flux conversion factors in the 0.5-2 
keV and in the 2-10 keV bands were computed using the re- 
sponse matrices at the aimpoint. A count rate of 1 count s _1 
in the soft band corresponds to a flux of (4.6 ± 0.1) x 10~ 12 
erg s _1 cm~ 2 and of (2.9 ± 0.3) x 10~ n erg s _1 cm~ 2 in the 
hard band. These numbers were computed assuming a Galactic 
absorbing column of 8 x 10 19 cm~ 2 and a power law spec- 
trum with a photon index T = 1.4 (i.e. the slope of the hard 
XRB). The quoted uncertainties are for photon indices in the 
range T = 1.1 — 1.7. As shown by Tozzi et al. (2001) a conver- 
sion factor changing with flux, in order to follow the hardening 
of the sources, results in a small variation of the Log N-Log S 
distribution, at a level of ~ 5%. 

A critical parameter in the search of faint sources is repre- 
sented by the detection threshold. In the context of detection 
by wavelet algorithms this is usually fixed in terms of spurious 
sources per field (Lazzati et al. 1998). To compare our results 
with those reported by Tozzi et al. (2001), we require to have 
less than 4.3 spurious sources within the selected area (i.e. 6 
spurious sources per Chandra field). This threshold has been 
accurately tested with simulations and the number of spurious 
sources is well under control (Moretti et al. 2001). This thresh- 
old corresponds to a source significance larger than 4 a. 

3. LOG N-LOG S RELATIONS IN THE SOFT AND HARD BANDS 

We searched the inner image with the BMW-Chandra de- 
tection algorithm described in Moretti et al. (2001). We de- 
tect 239 and 147 sources in the soft and hard band, respec- 
tively. There are 22 sources (~ 8% of the total number of de- 
tected sources) that are revealed only in the hard band, and 116 
sources (~ 44%) that are revealed only in the soft band. 

The sky coverage at a given flux, i.e. the area over which a 
given source can be detected above our threshold (as due to the 
varying vignetting, angular resolution and exposure time), has 
been derived by extensive simulations (100 fields per band) in 4 
circular annuli (see Moretti et al. 2001). The inhomogeneities 
of the ACIS-I field of view, such as gaps between the CCDs, 
are dealt with throrough simulations. 

The flux distribution of the detected sources suffers from the 
well known Eddington bias (e.g. Hasinger et al. 1993). Faint 
sources can be detected only if superposed on positive back- 
ground fluctuations (due either to a faint contaminating source 
or a statistical background fluctuation) and therefore their fluxes 
are systematically overestimated by factor as large as 2-3. This 
bias starts affecting the data at a level of ~ 20 counts in the 



soft band (~ 10 6 ergs -1 cm~ 2 ) and ~ 30 counts in the hard 
band (~ 10~ 15 ergs -1 cm -2 ). This effect has a stronger im- 
pact on the 2-10 keV Log N-Log S, due to the higher value 
of the background. Following the approach by Vikhlinin et al. 
(1995), we statistically correct for this bias. The completeness 
of the detection is then accounted for by the completeness func- 
tion. The accuracy of the procedure is assured by extensive 
simulations (see Moretti et al. 2001 and Figs. 4-5 therein). 
Our simulations show that we are able to recover the number 



source distribution down to 4 corrected counts (11 before the 
bias correction) in the inner core of the image, declining to 8 
corrected counts (14 before correction) in the outskirts for the 
soft band. For the hard band the limit on the corrected counts 
are 5 in the center and 10 in the outskirts. These counts give a 
flux limit in the inner region of 2 x 10~ 17 erg s _1 cm~ 2 in the 
soft band and 2 x 10~ 16 erg s _1 cm~ 2 in the hard band. 
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FIG. 1. — Soft band Log N-Log S from inner 8' radius of the 1 Ms ob- 
servation of the CDFS (thick solid line) in the range 2xl0~ 17 -3x 
10~ 14 erg cm~ 2 s _1 for an average V = 1.4 power-law spectrum. The 
gray thick line shows the la confidence region taking into account the num- 
ber statistics and the flux conversion uncertainties. The thick dashed solid line 
shows the best M-L fit power law, while the other lines show fits from previ- 
ous works. The thin solid line shows the ROSAT-HRI results on the Lockman 
hole (Hasinger et al. 1998); the shaded area shows the result of the background 
fluctuation analysis of ROSAT-HRI data (Hasinger et al. 1993), while the dot- 
ted and dash-dotted lines show the best fit from the analysis of the SSA13 
field Mushotzky et al. (2000) and of the first 120 ks on the CDFS Giacconi et 
al. (2001), respectively. The results on the first 300 ks of the CDFS (Tozzi et 
al. 2001) are not shown, since they would have been virtually indistinguish- 
able from our best fit line. The inset shows the 1, 2 and 3<r contours from 
the simultaneous M-L fit of the slope and normalization of the Log N-Log S. 
The asterisk, circle and diamond are the results of Mushotzky et al. (2000), 
Giacconi et al. (2001) and Tozzi et al. (2001), respectively. 



To characterize the flux distribution of sources detected in the 
soft band we performed a maximum likelihood fit to the differ- 
ential (unbinned) source flux distribution with a power-law. We 
obtain as a best fit: 



N(> S) = 360 



( S ) 

V2 x 10- 15 / 



-0.68 



CgS. 



(1) 



Fitting the same data, adaptively binning to 10 sources per bin, 
we obtain a good fit by means of a x 2 test (null hypothesis 
probability ~ 10%). This assures us of the goodness of the 
fit. The error at 90% confidence level for the normalization is 
Kis = 360lgy (if 15 means that fluxes are measured in units of 
2 x 10~ 15 erg s _1 cm -2 ), while the best fit value for the slope 
is a s = —0.68 ± 0.03 (see also Fig. 1). Both our normaliza- 
tion and slope, are in good agreement with ROSAT (Hasinger 
et al. 1998) Chandra (Tozzi et al. 2001, Rosati et al. 2001) and 
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XMM-Newton (Hasinger et al. 2001) up to their faintest flux 
levels. 

Due to Galactic contributions and following all previous 
studies, we consider the soft XRB in the 1-2 keV band only. 
Our fitted Log N-Log S distribution gives an integrated flux 
of ~ 5.4 x 10~ 13 erg s _1 cm~ 2 deg~ 2 for fluxes fainter than 
10~ 15 erg s _1 cm~ 2 . Summing the contribution at higher 
fluxes (Hasinger et al. 1998), the total flux of resolved point 
sources is <~ 3.5 x 10~ 12 erg s _1 cm~ 2 deg~ 2 . This accounts 
for 80% of the unresolved extragalactic background measured 
by ROSAT and by a joined ASCA/ROSAT analysis (Chen et al. 
1997). Extending to even lower fluxes the Log N-Log S dis- 
tribution (assuming the same slope) we can make up to 85% of 
the soft XRB, at the most. This is in line with recent predictions 
for the soft XRB produced by the warm/hot diffuse intergalac- 
tic medium which may contribute at a 5 — 15% level (Phillips et 
al. 2001). If, instead, we assume the lower (by a factor ~ 0.85) 
value measured by ASCA (Gendreau et al. 1995), then we are 
able to account for 96% of the flux in the soft XRB and fully 
account for it extending the source distribution at the lowest 
level. 

For the hard band our faint limit can be confidently taken 
down to 2 x 10~ 16 erg s _1 cm~ 2 (see Fig. 2). In order to 
be consistent with the results at larger fluxes, we normalized 
our distribution to the ASCA data (della Ceca et al. 2000) 3 . 
The Log N-Log S distribution cannot be described with a sin- 
gle power-law (probability of the null hypothesis <~ 0.02%, 
Fig. 2). We tried to model the data with a broken power-law 
and with a smoothly-jointed broken power-law. The second 
function gives a much better agreement with the data and it has 
the form: 



were the case a contribution from truly diffuse emission or from 
a new class of sources has to be advocated. 
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FIG. 2. — Hard band Log N-Log S from inner 8' radius of the 1 Ms 
observation of the CDFS (thick solid line) in the range 2 X 10~ 16 - 2 X 
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erg cm i s 1 for an average Y = 1.4 power-law spectrum. The dis- 



N(> S) = 1.2 x 10 4 



(2 x 10~ 15 ) Qhl 



cgs. (2) 



To match the bright flux end, we fixed the first slope to the 
ASCA value, i.e. a hl = 1.67 (della Ceca et al. 2000). 
The maximum likelihood fit yields a faint-end slope a n 2 = 



tribution is normalized to the ASCA value and slope (della Ceca et al. 2000). 
The gray thick line shows the lcr confidence region taking into account the 
number statistics and the flux conversion uncertainties. The thick dashed solid 
line shows the best M-L fit smoothly joined broken power law, while the 
other lines shows fits from previous works. The thin dashed line shows the 
BeppoSAX result (Vecchi et al. 1999) while the shaded area shows the result 
of the background fluctuation analysis of BeppoSAX data (Perri & Giommi 
2000). Lines related to Chandra fields are as in Fig. 1. The inset shows the 
1 , 2 and 3 cr contours from the simultaneous M-L fit of the faint slope a^2 
and break flux So , holding the bright slope fixed to the ASCA result and the 
normalization fixed to the best-fit value (see text). 



0.58 ± 0.03 and a break flux S = 1.7 



+0.3 
-0.2 



10 



-14 



erg s" 



cm 2 (all at 90% confidence level for a single parameter; see 
Fig. 2). 

The total contribution to the XRB is then 1.6 x 10~ n erg 
s _1 cm~ 2 deg~ 2 . This value has to be compared with the total 
(unresolved) XRB of 1.6 x 10~ n erg s _1 cm~ 2 deg~ 2 from 
UHURU and HEAO-1 (Marshall et al. 1980), as well as with 
the more recent determinations from the BeppoSAX and ASCA 
of 2.4 x 10~ n erg s _1 cm~ 2 deg~ 2 (Chiappetti et al. 1998; 
Vecchi et al. 1999; Ueda et al. 1999; Ishisaki et al. 2000; Perri 
& Giommi 2000), in agreement with the old Wisconsin mea- 
surements (McCammon et al. 1983). Our results are in good 
agreement with the UHURU/HEAO- 1 value (considering also 
the error in the flux determination due to the uncertainties in 
the background slope, which amounts to 10%). The slope of 
the faint end counts converges slower than logarithmically and 
thus the faint sources provide a small contribution to the X-ray 
background (adding only 10% more extrapolating down to zero 
flux). In this case we can claim to have fully resolved the 2-10 
keV hard XRB. On the contrary, if we adopt the BeppoSAX and 
ASCA value, the fraction of resolved sources in the CDFS low- 
ers to 68%. If we extrapolate our best fit distribution to lower 
fluxes, we can make up 73% of the entire XRB, at most. If this 

3 This normalization of the Log N-Log S distribution allows to get rid of cosmic variance at bright fluxes. At faint fluxes this is not true and the small field analyzed 
can play likely a role, e.g. in the Chandra observation of the HDF more sources are detected at the same flux level (Brandt et al. 2001). 



4. CONCLUSIONS 

The origin of the XRB has been extensively studied since its 
discovery (Giacconi et al. 1962). In the case of the soft (1-2 
keV) XRB there is general consensus that the bulk of the emis- 
sion can be ascribed to broad line AGN (i.e. Seyfert 1 galaxies; 
Hasinger et al. 1998; Schmidt et al. 1998). At faint fluxes 
( ^ 10~ 15 erg s _1 cm~ 2 ) nearby (z ^ 0.6) optically-normal 
galaxies are also being detected as soft sources (Barger et al. 
2001; Tozzi et al. 2001; Brandt et al. 2001; Schreier et al. 
2001). Our analysis of the CDFS extends these studies down 
to 2 x 10 -17 erg s _1 cm~ 2 . We found that even at these low 
fluxes the Log N-Log S distribution can still be represented by 
the extrapolation from higher fluxes (Tozzi et al. 2001) and we 
are able to resolve in point sources > 80% of the soft XRB, the 
main uncertainties being the sky flux itself. 

In the hard band the XRB is made by the superposition of 
absorbed and unabsorbed sources (Setti & Woltjer 1989). In 
particular, the flat spectrum of the background (V ~ 1.4) im- 
plies a considerable absorption in most objects. Several studies 
have emphasized the importance of the hard XRB as a probe 
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of the evolution of the black hole population (Fabian & Iwa- 
sawa 1999). The small number of hard sources not detected in 
the soft band indicate that even in highly obscured AGN a size- 
able soft emission can still be produced due to, e.g., scattering, 
partial covering of the central radiation or from starburst emis- 
sion around the AGN (e.g. Turner et al. 1997). This confirms 
that also at these extreme flux levels the ratio between soft and 
hard X-ray emission is ~ 1 — 10%. In addition, as shown in 
Barger et al. (2001), the effects of obscuration get significantly 
reduced at z ~ 2, since we observe in the soft band the harder 
1 .5-6 keV rest frame interval. 

When matching the Log N-Log S distribution from 
ASCA/BeppoSAX with deeper data the presence of a knee is 
evident. We are able to constrain the Log N-Log S relation 
below and above the knee, which occurs at ~ 2 x 10~ 14 erg 
s _1 cm~ 2 . We found a fit to the data with a smoothly joined 
power law with the bright end fixed to the ASCA value (1.67) 
and the fainter one to ah2 ~ 0.6. At a level of ~ 2 x 10~ 16 erg 



s _1 cm~ 2 we are fully consistent with the UHURU/HEAO-1 
estimate, whereas we make up ~ 68% of the hard XRB, if we 
assume the BeppoSAX/ASCA value. 

A dedicated, multiwavelenght follow up is mandatory to 
exploit the enormous potential of this (and similar) datasets 
(Barger et al. 2001; Giacconi et al. 2001b). The identifica- 
tion of X-ray sources at the faintest limit and their optical char- 
acterization (by means of photometric redshift) can solve the 
issue of the evolution of obscured AGN population at very high 
redshift. 
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